A new antibiotic traps lipopolysaccharide in its intermembrane transporter.

Gram-negative bacteria are extraordinarily difficult to kill because their cytoplasmic membrane is surrounded by an outer membrane that blocks the entry of most antibiotics. The impenetrable nature of the outer membrane is due to the presence of a large, amphipathic glycolipid called lipopolysaccharide (LPS) in its outer leaflet1. Assembly of the outer membrane requires transport of LPS across a protein bridge that spans from the cytoplasmic membrane to the cell surface. Maintaining outer membrane integrity is essential for bacterial cell viability, and its disruption can increase susceptibility to other antibiotics2-6. Thus, inhibitors of the seven lipopolysaccharide transport (Lpt) proteins that form this transenvelope transporter have long been sought. A new class of antibiotics that targets the LPS transport machine in Acinetobacter was recently identified. Here, using structural, biochemical and genetic approaches, we show that these antibiotics trap a substrate-bound conformation of the LPS transporter that stalls this machine. The inhibitors accomplish this by recognizing a composite binding site made up of both the Lpt transporter and its LPS substrate. Collectively, our findings identify an unusual mechanism of lipid transport inhibition, reveal a druggable conformation of the Lpt transporter and provide the foundation for extending this class of antibiotics to other Gram-negative pathogens.

V3 tip determinants of susceptibility to inhibition by CD4-mimetic compounds in natural clade A human immunodeficiency virus (HIV-1) envelope glycoproteins.

IMPORTANCE: CD4-mimetic compounds (CD4mcs) are small-molecule inhibitors of human immunodeficiency virus (HIV-1) entry into host cells. CD4mcs target a pocket on the viral envelope glycoprotein (Env) spike that is used for binding to the receptor, CD4, and is highly conserved among HIV-1 strains. Nonetheless, naturally occurring HIV-1 strains exhibit a wide range of sensitivities to CD4mcs. Our study identifies changes distant from the binding pocket that can influence the susceptibility of natural HIV-1 strains to the antiviral effects of multiple CD4mcs. We relate the antiviral potency of the CD4mc against this panel of HIV-1 variants to the ability of the CD4mc to activate entry-related changes in Env conformation prematurely. These findings will guide efforts to improve the potency and breadth of CD4mcs against natural HIV-1 variants.

Cryo-EM structures reveal native GABAA receptor assemblies and pharmacology.

Type A γ-aminobutyric acid receptors (GABAARs) are the principal inhibitory receptors in the brain and the target of a wide range of clinical agents, including anaesthetics, sedatives, hypnotics and antidepressants1-3. However, our understanding of GABAAR pharmacology has been hindered by the vast number of pentameric assemblies that can be derived from 19 different subunits4 and the lack of structural knowledge of clinically relevant receptors. Here, we isolate native murine GABAAR assemblies containing the widely expressed α1 subunit and elucidate their structures in complex with drugs used to treat insomnia (zolpidem (ZOL) and flurazepam) and postpartum depression (the neurosteroid allopregnanolone (APG)). Using cryo-electron microscopy (cryo-EM) analysis and single-molecule photobleaching experiments, we uncover three major structural populations in the brain: the canonical α1ß2γ2 receptor containing two α1 subunits, and two assemblies containing one α1 and either an α2 or α3 subunit, in which the single α1-containing receptors feature a more compact arrangement between the transmembrane and extracellular domains. Interestingly, APG is bound at the transmembrane α/ß subunit interface, even when not added to the sample, revealing an important role for endogenous neurosteroids in modulating native GABAARs. Together with structurally engaged lipids, neurosteroids produce global conformational changes throughout the receptor that modify the ion channel pore and the binding sites for GABA and insomnia medications. Our data reveal the major α1-containing GABAAR assemblies, bound with endogenous neurosteroid, thus defining a structural landscape from which subtype-specific drugs can be developed.

Molecular mechanisms of SARS-CoV-2 resistance to nirmatrelvir.

Nirmatrelvir is a specific antiviral drug that targets the main protease (Mpro) of SARS-CoV-2 and has been approved to treat COVID-191,2. As an RNA virus characterized by high mutation rates, whether SARS-CoV-2 will develop resistance to nirmatrelvir is a question of concern. Our previous studies have shown that several mutational pathways confer resistance to nirmatrelvir, but some result in a loss of viral replicative fitness, which is then compensated for by additional alterations3. The molecular mechanisms for this observed resistance are unknown. Here we combined biochemical and structural methods to demonstrate that alterations at the substrate-binding pocket of Mpro can allow SARS-CoV-2 to develop resistance to nirmatrelvir in two distinct ways. Comprehensive studies of the structures of 14 Mpro mutants in complex with drugs or substrate revealed that alterations at the S1 and S4 subsites substantially decreased the level of inhibitor binding, whereas alterations at the S2 and S4' subsites unexpectedly increased protease activity. Both mechanisms contributed to nirmatrelvir resistance, with the latter compensating for the loss in enzymatic activity of the former, which in turn accounted for the restoration of viral replicative fitness, as observed previously3. Such a profile was also observed for ensitrelvir, another clinically relevant Mpro inhibitor. These results shed light on the mechanisms by which SARS-CoV-2 evolves to develop resistance to the current generation of protease inhibitors and provide the basis for the design of next-generation Mpro inhibitors.

High-resolution landscape of an antibiotic binding site.

Antibiotic binding sites are located in important domains of essential enzymes and have been extensively studied in the context of resistance mutations; however, their study is limited by positive selection. Using multiplex genome engineering1 to overcome this constraint, we generate and characterize a collection of 760 single-residue mutants encompassing the entire rifampicin binding site of Escherichia coli RNA polymerase (RNAP). By genetically mapping drug-enzyme interactions, we identify an alpha helix where mutations considerably enhance or disrupt rifampicin binding. We find mutations in this region that prolong antibiotic binding, converting rifampicin from a bacteriostatic to bactericidal drug by inducing lethal DNA breaks. The latter are replication dependent, indicating that rifampicin kills by causing detrimental transcription-replication conflicts at promoters. We also identify additional binding site mutations that greatly increase the speed of RNAP.Fast RNAP depletes the cell of nucleotides, alters cell sensitivity to different antibiotics and provides a cold growth advantage. Finally, by mapping natural rpoB sequence diversity, we discover that functional rifampicin binding site mutations that alter RNAP properties or confer drug resistance occur frequently in nature.

MEN1 mutations mediate clinical resistance to menin inhibition.

Chromatin-binding proteins are critical regulators of cell state in haematopoiesis1,2. Acute leukaemias driven by rearrangement of the mixed lineage leukaemia 1 gene (KMT2Ar) or mutation of the nucleophosmin gene (NPM1) require the chromatin adapter protein menin, encoded by the MEN1 gene, to sustain aberrant leukaemogenic gene expression programs3-5. In a phase 1 first-in-human clinical trial, the menin inhibitor revumenib, which is designed to disrupt the menin-MLL1 interaction, induced clinical responses in patients with leukaemia with KMT2Ar or mutated NPM1 (ref. 6). Here we identified somatic mutations in MEN1 at the revumenib-menin interface in patients with acquired resistance to menin inhibition. Consistent with the genetic data in patients, inhibitor-menin interface mutations represent a conserved mechanism of therapeutic resistance in xenograft models and in an unbiased base-editor screen. These mutants attenuate drug-target binding by generating structural perturbations that impact small-molecule binding but not the interaction with the natural ligand MLL1, and prevent inhibitor-induced eviction of menin and MLL1 from chromatin. To our knowledge, this study is the first to demonstrate that a chromatin-targeting therapeutic drug exerts sufficient selection pressure in patients to drive the evolution of escape mutants that lead to sustained chromatin occupancy, suggesting a common mechanism of therapeutic resistance.

Direct Blockade of the Norovirus Histo-Blood Group Antigen Binding Pocket by Nanobodies.

Noroviruses are the leading cause of outbreaks of acute gastroenteritis. These viruses usually interact with histo-blood group antigens (HBGAs), which are considered essential cofactors for norovirus infection. This study structurally characterizes nanobodies developed against the clinically important GII.4 and GII.17 noroviruses with a focus on the identification of novel nanobodies that efficiently block the HBGA binding site. Using X-ray crystallography, we have characterized nine different nanobodies that bound to the top, side, or bottom of the P domain. The eight nanobodies that bound to the top or side of the P domain were mainly genotype specific, while one nanobody that bound to the bottom cross-reacted against several genotypes and showed HBGA blocking potential. The four nanobodies that bound to the top of the P domain also inhibited HBGA binding, and structural analysis revealed that these nanobodies interacted with several GII.4 and GII.17 P domain residues that commonly engaged HBGAs. Moreover, these nanobody complementarity-determining regions (CDRs) extended completely into the cofactor pockets and would likely impede HBGA engagement. The atomic level information for these nanobodies and their corresponding binding sites provide a valuable template for the discovery of additional "designer" nanobodies. These next-generation nanobodies would be designed to target other important genotypes and variants, while maintaining cofactor interference. Finally, our results clearly demonstrate for the first time that nanobodies directly targeting the HBGA binding site can function as potent norovirus inhibitors. IMPORTANCE Human noroviruses are highly contagious and a major problem in closed institutions, such as schools, hospitals, and cruise ships. Reducing norovirus infections is challenging on multiple levels and includes the frequent emergence of antigenic variants, which complicates designing effective, broadly reactive capsid therapeutics. We successfully developed and characterized four norovirus nanobodies that bound at the HBGA pockets. Compared with previously developed norovirus nanobodies that inhibited HBGA through disrupted particle stability, these four novel nanobodies directly inhibited HBGA engagement and interacted with HBGA binding residues. Importantly, these new nanobodies specifically target two genotypes that have caused the majority of outbreaks worldwide and consequently would have an enormous benefit if they could be further developed as norovirus therapeutics. To date, we have structurally characterized 16 different GII nanobody complexes, a number of which block HBGA binding. These structural data could be used to design multivalent nanobody constructs with improved inhibition properties.

Probing the Requirements for Dual Angiotensin-Converting Enzyme C-Domain Selective/Neprilysin Inhibition.

Selective inhibition of the angiotensin-converting enzyme C-domain (cACE) and neprilysin (NEP), leaving the ACE N-domain (nACE) free to degrade bradykinin and other peptides, has the potential to provide the potent antihypertensive and cardioprotective benefits observed for nonselective dual ACE/NEP inhibitors, such as omapatrilat, without the increased risk of adverse effects. We have synthesized three 1-carboxy-3-phenylpropyl dipeptide inhibitors with nanomolar potency based on the previously reported C-domain selective ACE inhibitor lisinopril-tryptophan (LisW) to probe the structural requirements for potent dual cACE/NEP inhibition. Here we report the synthesis, enzyme kinetic data, and high-resolution crystal structures of these inhibitors bound to nACE and cACE, providing valuable insight into the factors driving potency and selectivity. Overall, these results highlight the importance of the interplay between the S1' and S2' subsites for ACE domain selectivity, providing guidance for future chemistry efforts toward the development of dual cACE/NEP inhibitors.

Interaction between miR4749 and Human Serum Albumin as Revealed by Fluorescence, FRET, Atomic Force Spectroscopy and Computational Modelling.

The interaction of Human Serum Albumin (HSA) with the microRNA, miR4749, was investigated by Atomic Force Spectrscopy (AFS), static and time-resolved fluorescence spectroscopy and by computational methods. The formation of a HSA/miR4749 complex with an affinity of about 104 M-1 has been assessed through a Stern-Volmer analysis of steady-state fluorescence quenching of the lone Trp residue (Trp214) emission of HSA. Förster Resonance Energy Transfer (FRET) measurements of fluorescence lifetime of the HSA/miR4749 complex were carried out in the absence and in the presence of an acceptor chromophore linked to miR4749. This allowed us to determine a distance of 4.3 ± 0.5 nm between the lone Trp of HSA and the dye bound to miR4749 5p-end. Such a distance was exploited for a screening of the possible binding sites between HSA and miR4749, as predicted by computational docking. Such an approach, further refined by binding free energy calculations, led us to the identification of a consistent model for the structure of the HSA/miR4749 complex in which a positively charged HSA pocket accommodates the negatively charged miRNA molecule. These results designate native HSA as a suitable miRNA carrier under physiological conditions for delivering to appropriate targets.

In Silico and Experimental ADAM17 Kinetic Modeling as Basis for Future Screening System for Modulators.

Understanding the mechanisms of modulators' action on enzymes is crucial for optimizing and designing pharmaceutical substances. The acute inflammatory response, in particular, is regulated mainly by a disintegrin and metalloproteinase (ADAM) 17. ADAM17 processes several disease mediators such as TNFα and APP, releasing their soluble ectodomains (shedding). A malfunction of this process leads to a disturbed inflammatory response. Chemical protease inhibitors such as TAPI-1 were used in the past to inhibit ADAM17 proteolytic activity. However, due to ADAM17's broad expression and activity profile, the development of active-site-directed ADAM17 inhibitor was discontinued. New 'exosite' (secondary substrate binding site) inhibitors with substrate selectivity raised the hope of a substrate-selective modulation as a promising approach for inflammatory disease therapy. This work aimed to develop a high-throughput screen for potential ADAM17 modulators as therapeutic drugs. By combining experimental and in silico methods (structural modeling and docking), we modeled the kinetics of ADAM17 inhibitor. The results explain ADAM17 inhibition mechanisms and give a methodology for studying selective inhibition towards the design of pharmaceutical substances with higher selectivity.

Angiogenesis and anti-leukaemia activity of novel indole derivatives as potent colchicine binding site inhibitors.

The screened compound DYT-1 from our in-house library was taken as a lead (inhibiting tubulin polymerisation: IC50=25.6 µM, anti-angiogenesis in Zebrafish: IC50=38.4 µM, anti-proliferation against K562 and Jurkat: IC50=6.2 and 7.9 µM, respectively). Further investigation of medicinal chemistry conditions yielded compound 29e (inhibiting tubulin polymerisation: IC50=4.8 µM and anti-angiogenesis in Zebrafish: IC50=3.6 µM) based on tubulin and zebrafish assays, which displayed noteworthily nanomolar potency against a variety of leukaemia cell lines (IC50= 0.09-1.22 µM), especially K562 cells where apoptosis was induced. Molecular docking, molecular dynamics (MD) simulation, radioligand binding assay and cellular microtubule networks disruption results showed that 29e stably binds to the tubulin colchicine site. 29e significantly inhibited HUVEC tube formation, migration and invasion in vitro. Anti-angiogenesis in vivo was confirmed by zebrafish xenograft. 29e also prominently blocked K562 cell proliferation and metastasis in blood vessels and surrounding tissues of the zebrafish xenograft model. Together with promising physicochemical property and metabolic stability, 29e could be considered an effective anti-angiogenesis and -leukaemia drug candidate that binds to the tubulin colchicine site.

Discovery and Characterization of a Cryptic Secondary Binding Site in the Molecular Chaperone HSP70.

Heat Shock Protein 70s (HSP70s) are key molecular chaperones that are overexpressed in many cancers and often associated with metastasis and poor prognosis. It has proven difficult to develop ATP-competitive, drug-like small molecule inhibitors of HSP70s due to the flexible and hydrophilic nature of the HSP70 ATP-binding site and its high affinity for endogenous nucleotides. The aim of this study was to explore the potential for the inhibition of HSP70 through alternative binding sites using fragment-based approaches. A surface plasmon resonance (SPR) fragment screen designed to detect secondary binding sites in HSP70 led to the identification by X-ray crystallography of a cryptic binding site in the nucleotide-binding domain (NBD) of HSP70 adjacent to the ATP-binding site. Fragment binding was confirmed and characterized as ATP-competitive using SPR and ligand-observed NMR methods. Molecular dynamics simulations were applied to understand the interactions with the protein upon ligand binding, and local secondary structure changes consistent with interconversion between the observed crystal structures with and without the cryptic pocket were detected. A virtual high-throughput screen (vHTS) against the cryptic pocket was conducted, and five compounds with diverse chemical scaffolds were confirmed to bind to HSP70 with micromolar affinity by SPR. These results identified and characterized a new targetable site on HSP70. While targeting HSP70 remains challenging, the new site may provide opportunities to develop allosteric ATP-competitive inhibitors with differentiated physicochemical properties from current series.

Metal Ion Binding Induces Local Protein Unfolding and Destabilizes Human Carbonic Anhydrase II.

Human carbonic anhydrase II (HCA) is a robust metalloprotein and an excellent biological model system to study the thermodynamics of metal ion coordination. Apo-HCA binds one zinc ion or two copper ions. We studied these binding processes at five temperatures (15-35 °C) using isothermal titration calorimetry, yielding thermodynamic parameters corrected for pH and buffer effects. We then sought to identify binding-induced structural changes. Our data suggest that binding at the active site organizes 6-8 residues; however, copper binding near the N-terminus results in a net unfolding of 6-7 residues. This surprising destabilization was confirmed using circular dichroism and protein stability measurements. Metal binding induced unfolding may represent an important regulatory mechanism, but it may be easily missed by NMR and X-ray crystallography. Thus, in addition to highlighting a highly novel binding-induced unfolding event, we demonstrate the value of calorimetry for studying the structural implications of metal binding.

Identification of novel non-toxic and anti-angiogenic α-fluorinated chalcones as potent colchicine binding site inhibitors.

α-Fluorinated chalcones were prepared and evaluated for their cell growth inhibitory properties against six human cancer cell lines. The most potent chalcone 4c demonstrated excellent selective toxicity against cancer cells versus normal human cells, with IC50 values at nanomolar concentration ranges against 5 cancer cell lines. A further study revealed that 4c could bind to the colchicine site of tubulin, disrupt the cell microtubule networks, and effectively inhibit tubulin polymerisation. Cellular-based mechanism studies elucidated that 4c arrested MGC-803 cell cycle at G2/M phase. In addition, 4c dose-dependently caused Caspase-induced apoptosis of MGC-803 cells through mitochondrial dysfunction. Notably, compound 4c was found to inhibit the HUVECs tube formation, migration, and invasion in vitro. Furthermore, our data suggested that treatment with 4c significantly reduced MGC-803 cells metastasis and proliferation in vitro. Overall, this work showed that chalcone hybrid 4c is a potent inhibitor of tubulin assembly with prominent anti-angiogenesis and anti-cancer properties.

Cysteine-Based Protein Covalent Binding and Hepatotoxicity Induced by Emodin.

Emodin (EMD) is a major ingredient of Polygonum multiflorum Thunb. (PMT), which has shown adverse liver reactions. Despite multiple pharmacological activities, EMD is reported to show various toxicities. Our early study demonstrated the reactivity of EMD to glutathione. This study aimed to determine the covalent interaction of hepatic protein with EMD and the correlation of the protein modification with hepatotoxicity induced by EMD. EMD-derived protein adduction was detected in an incubation mixture containing mouse liver homogenates and EMD. Such protein adduction was also observed in hepatic protein obtained from mice exposed to EMD. The protein covalent binding occurred in time- and dose-dependent manners. Pre-treatment of l-buthionine-sulfoximine significantly potentiated EMD-induced adduction and hepatotoxicity caused by EMD and lipopolysaccharide co-treatment. As expected, EMD-derived protein modification was observed in mouse primary hepatocytes treated with EMD. The increase in EMD exposure concentration intensified EMD-derived protein adduction and increased EMD-induced cell death. The susceptibility of hepatocytes to EMD cytotoxicity and the intensity of EMD-induced protein adduction were attenuated by the co-treatment of hepatocytes with N-acetyl cysteine. A good association of protein modification with hepatotoxicity induced by EMD was illustrated, which facilitates the understanding of the mechanism of hepatotoxicity induced by EMD.

Amyloidogenic immunoglobulin light chain kinetic stabilizers comprising a simple urea linker module reveal a novel binding sub-site.

In immunoglobulin light chain (LC) amyloidosis, the misfolding, or misfolding and misassembly of LC a protein or fragments thereof resulting from aberrant endoproteolysis, causes organ damage to patients. A small molecule "kinetic stabilizer" drug could slow or stop these processes and improve prognosis. We previously identified coumarin-based kinetic stabilizers of LCs that can be divided into four components, including a "linker module" and "distal substructure". Our prior studies focused on characterizing carbamate, hydantoin, and spirocyclic urea linker modules, which bind in a solvent-exposed site at the VL-VL domain interface of the LC dimer. Here, we report structure-activity relationship data on 7-diethylamino coumarin-based kinetic stabilizers. This substructure occupies the previously characterized "anchor cavity" and the "aromatic slit". The potencies of amide and urea linker modules terminating in a variety of distal substructures attached at the 3-position of this coumarin ring were assessed. Surprisingly, crystallographic data on a 7-diethylamino coumarin-based kinetic stabilizer reveals that the urea linker module and distal substructure attached at the 3-position bind a solvent-exposed region of the full-length LC dimer distinct from previously characterized sites. Our results further elaborate the small-molecule binding surface of LCs that could be occupied by potent and selective LC kinetic stabilizers.

Alteration of a Cry1A Shared Binding Site in a Cry1Ab-Selected Colony of Ostrinia furnacalis.

The Asian corn borer, Ostrinia furnacalis (Guenée, 1854), is a highly damaging pest in Asia and the Pacific islands, and larvae feed mainly from corn crops. To determine the suitability of Bt-corn technology for the future control of this pest, understanding the potential to develop resistance to Cry1Ab and the basis of cross-resistance to other Cry1 proteins is of great interest. Here, we have explored the binding of Cry1A proteins to brush border membrane vesicles from two O. furnacalis colonies, one susceptible (ACB-BtS) and one laboratory-selected with Cry1Ab (ACB-AbR). The insects developed resistance to Cry1Ab and showed cross-resistance to Cry1Aa, Cry1Ac, and Cry1F. Binding assays with radiolabeled Cry1Ab and brush border membrane vesicles from susceptible insects showed that Cry1A proteins shared binding sites, though the results were not conclusive for Cry1F. The results were confirmed using radiolabeled Cry1Aa. The resistant insects showed a reduction of the specific binding of both Cry1Ab and Cry1Aa, suggesting that part of the binding sites were lost or altered. Competition binding assays showed full competition between Cry1Ab and Cry1Aa proteins in the susceptible colony but only partial competition in resistant insects, confirming the alteration of some, but not all, binding sites for these two proteins. The binding site model for Cry1A proteins in O. furnacalis is in agreement with the occurrence of multiple membrane receptors for these proteins.

Discovery of Pyrazolo[1,5-a]pyrazin-4-ones as Potent and Brain Penetrant GluN2A-Selective Positive Allosteric Modulators Reducing AMPA Receptor Binding Activity.

N-Methyl-d-aspartate receptors (NMDARs) are members of the ionotropic glutamate receptor family and play a crucial role in learning and memory by regulating synaptic plasticity. Activation of NMDARs containing GluN2A, one of the NMDAR subunits, has recently attracted attention as a promising therapeutic approach for neuropsychiatric diseases such as schizophrenia, depression, and epilepsy. In the present study, we developed potent and brain-penetrable GluN2A-selective positive allosteric modulators. Lead compound 2b was generated by scaffold hopping of hit compound 1, identified from the internal alpha-amino-3-hydroxy-5-methyl-4-isoxazolepropionic acid receptor (AMPAR)-focused compound library through a high-throughput screening campaign. Subsequent optimization of the lead compound, including a structure-based drug design approach, resulted in the identification of a potent GluN2A PAM (R)-9, which possessed high selectivity against both subtypes of AMPAR and NMDAR. Furthermore, (R)-9 significantly enhanced long-term potentiation in the rat hippocampus 24 h after oral administration, indicating that this molecule is a potentially useful in vivo pharmacological tool for treating psychiatric diseases.

Advances in Magnetic Microbead Affinity Selection Screening: Discovery of Natural Ligands to the SARS-CoV-2 Spike Protein.

Affinity selection-mass spectrometry, which includes magnetic microbead affinity selection-screening (MagMASS), is ideal for the discovery of ligands in complex mixtures that bind to pharmacological targets. Therapeutic agents are needed to prevent or treat COVID-19, which is caused by the severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). Infection of human cells by SARS-CoV-2 involves binding of the virus spike protein subunit 1 (S1) to the human cell receptor angiotensin converting enzyme-2 (ACE2). Like antibodies, small molecules have the potential to block the interaction of the viral S1 protein with human ACE2 and prevent SARS-CoV-2 infection. Therefore, a MagMASS assay was developed for the discovery of ligands to the S1 protein. Unlike previous MagMASS approaches, this new assay used robotics for 5-fold enhancement of throughput and sensitivity. The assay was validated using the SBP-1 peptide, which is identical to the ACE2 amino acid sequence recognized by the S1 protein, and then applied to the discovery of natural ligands from botanical extracts. Small molecule ligands to the S1 protein were discovered in extracts of the licorice species, Glycyrrhiza inflata. In particular, the licorice ligand licochalcone A was identified through dereplication and comparison with standards using HPLC with high-resolution tandem mass spectrometry.

Reconstruction of the binding pathway of an anti-HIV drug, Indinavir, in complex with the HTLV-1 protease using unaggregated unbiased molecular dynamics simulation.

Retroviruses are a growing concern for the health of human beings, and one of the dangerous members of this family is the Human T-cell Leukemia Virus 1 (HTLV-1) virus. It has affected more than 20 million people so far, and since there are no registered treatments against it yet, urgent treatment solutions are needed. One of the most promising drug targets to fight this virus is the protease enzyme of the virus's protein machinery. In this study, by utilizing a computational method called Unaggregated Unbiased Molecular Dynamics (UUMD), we reconstructed the binding pathway of a HTLV-1 protease inhibitor, Indinavir, to find the details of the binding pathway, the influential residues, and also the stable states of the binding pathway. We achieved the native conformation of the inhibitor in 6 rounds, 360 replicas by performing over 4 micro-seconds of UMD simulations. We found 3 Intermediate states between the solvated state and the native conformation state in the binding pathway. We also discovered that aromatic residues such as Trp98 and Trp98', catalytic residues Asp32 and Asp32', and the flap region's residues have the most influential roles in the binding pathway and also have the most contribution to the total interaction energies. We believe that the details found in this study would be a great guide for developing new treatment solutions against the HTLV-1 virus by inhibiting the HTLV-1 protease.